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Abstract: We consider a class of eigenvalue problems for poly-harmonic oper- 
ators, including Dirichlet and buckling-type eigenvalue problems. We prove an 
analyticity result for the dependence of the symmetric functions of the eigenval- 
ues upon domain perturbations and compute Hadamard-type formulas for the 
Frechet differentials. We also consider isovolumetric domain perturbations and 
P_( ■ characterize the corresponding critical domains for the symmetric functions of 

, the eigenvalues. Finally, we prove that balls are critical domains. 
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tj- ' 1 Introduction 

Let ?i,mGNo with < m < n and O be a bounded open set in M. N with smooth 
^ | boundary. We consider the following eigenvalue problem 

(N 

-A) n u = A(-A) m u, in fi, 

0, on d£l, 




where v denotes the unit outer normal to dfl. The case m = corresponds to the 
well-known eigenvalue problem for the poly-harmonic operator (— A) n subject 
to Dirichlet boundary conditions, while the case m > represents a buckling- 
type problem. These cases include important problems in linear elasticity. For 
instance, for N = 2, Vio arises in the study of a vibrating membrane stretched in 
a fixed frame, V20 corresponds to the case of a vibrating clamped plate and V21 
is related to plate buckling. We are mainly interested in the Dirichlet problem 
V n o and the buckling problem P 21 . However, we prefer to present a unified 
approach involving all cases. Problem V nm admits a divergent sequence of positive 
eigenvalues of finite multiplicity represented as follows 

< Ai[fi] < X 2 M < ■ < Xjp] < ... . 

As usual, we agree to repeat each eigenvalue as many times as its multiplicity. 

In this paper we are interested in the dependence of Aj[fi] on Q. There is a 
vast literature devoted to domain perturbation problems for elliptic operators. 
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In particular, the cases n = 1, 2 which correspond to the Laplace operator and 
the bi-harmonic operator respectively, have been intensively studied by many 
authors. We refer to Bucur and Buttazzo |6j, Daners [9], Hale [12], Henry [13], 
Henrot [15], Kesavan [21] for updated information on this topic. The case n > 2 
has been much less investigated. However, a renewed general interest in higher 
order operators has been growing in the last decades as it appears in the extensive 
monograph by Gazzola, Grunau and Sweers [10] devoted to recent developments 
in the theory of poly-harmonic operators. As for domain perturbation problems, 
we refer to the papers [3] and [1] where spectral stability estimates for elliptic 
operators of arbitrary order are proved. See also the survey paper [5]. 

Our work is inspired by classical problems in spectral optimization, in par- 
ticular by the celebrated Rayleigh conjecture (see e.g., [IS]). Recall that the 
Rayleigh-Faber-Krahn inequality states that the first eigenvalue of the Laplace 
operator with Dirichlet boundary conditions (problem Vio) is minimized by the 
ball in the class of all bounded domains with a fixed measure. In symbols, 

Ai[fi*] < Aiffi], (1.2) 

where Q* is a ball with the same measure of Q. As for the bi-harmonic operator 
with Dirichlet boundary condition (problem V20), inequality (II. 2p was proved 
by Nadirashvili [20] for N = 2 and by Ashbaugh and Benguria [2J for N = 
3. See also Mohr [19]. Inequality (II. 2p can be proved also for plate buckling 
(problem P 2 i) under some extra assumptions, see [T5]. It should be noted that 
not much is known for higher eigenvalues Aj[f2] for j > 2. As a corollary of a 
general result by Buttazzo and Dal Maso [7], it is known that each eigenvalue 
of the Dirichlet Laplacian admits a minimizer in the class of all quasi-open sets 
with fixed measure, contained in a prescribed bounded region. However, that 
result says little about the shape of such minimizer. It is proved in Wolf and 
Keller [21] that the minimizers of higher eigenvalues in general are not balls and 
not even unions of balls. Moreover, by looking at the interesting numerical results 
presented in Oudet [23], one may get the idea that balls are not much relevant 
in the analysis of higher eigenvalues. 

Our main aim is to point out that, despite the above mentioned negative 
results, balls play a relevant role in the study of isovolumetric perturbations of Q 
for all eigenvalues Aj[fi] of all problems V nm . To do so we shall consider problems 
Vnm on families of open sets 0(f2) described by suitable diffeomorphisms defined 
on a fixed open set Q and we shall study the dependence of Aj[0(f2)] on 0. As 
is known, this allows to talk about differentiability and apply calculus in order 
to find critical eigenvalues with respect to perturbations of 0. One of the main 
difficulties in the analysis of higher eigenvalues is related to the variation of their 
multiplicity. This leads to bifurcation phenomena which complicate things. For 
instance, if a fixed open set 0(fi) is subject to a perturbation of then a 
multiple eigenvalue Aj[0(f2)] of multiplicity r may split into r simple eigenvalues 
Aj[0(f2)], . . . , Aj +r _i[0(f2)] in such a way that \j[<f)(Q)], . . . , Aj +r _i[0(f2)] are not 
different iable in the variable 0. As for the Laplace operator with Dirichlet or 
Neumann boundary conditions, it was pointed out in [16] that in the case of 
multiple eigenvalues it is natural to consider the elementary symmetric functions 
of the eigenvalues Aj[0(f2)], . . . , A J+r _i[0(f2)]. In this paper, we generalize the 
results of [HJ [T7] to all problems V nm . Namely, we prove that the elementary 
symmetric functions of the eigenvalues Aj[0(O)], . . . , Aj +r _i[0(fi)] of V nm depend 
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real analytically on (Theorem I3.2p and we prove that if 0(0) is a ball then 
is a critical point for such functions under volume constraint (Theorem 14. 7p . In 
fact, all critical points for the symmetric functions of the eigenvalues splitting 
from an eigenvalue A of multiplicity r, can be characterized as those open sets for 
which the following overdetermined system has a nontrivial solution (ui, . . . , u r ) 

{-A) n Ui = X(-A) m Ui, in 0(0), Vi = l,...,r, 

Ui = % = ..- = f^ = 0, ond0(O), Vi = l,...,r, (1.3) 
(f^) 2 + ••■ + (f^) 2 = const, on 50(0). 

Since (jl.3p is satisfied if 0(0) is a ball, it would be interesting to clarify whether 
the existence of nontrivial solutions to ( 11. 3ft on bounded connected open set 0(0) 
implies that 0(0) is a ball. For r = l, n = l,m = or n = 2,m = this can 
be proved by using the celebrated moving plane method under the assumption 
that the solution u\ does not change sign (see e.g., Henry [14] for the Laplace 
operator and Dalmasso [8] for the biharmonic operator); for r = 1, n = 2,m = 1 
a different method by Weinberger and Willms leads to the same conclusion (see 
e.g., [IB]). 



2 Notation and preliminaries 

Let n G N and O be a bounded open set in K. . By W n,2 (fl) we denote the 
Sobolev space of all functions in L 2 (0) which admit weak derivatives in L 2 (0) 
up to order n. By W^' 2 (Q) we denote the closure in W n,2 (Q) of the space of 
C^-functions with compact support in O. We consider the weak formulation of 
problem ( II. ip . To do so, for any m G No with < m < n, we consider the 
poly-harmonic operator A m as the operator from Wq ,2 {VL) to its dual (Wq ,2 (Q))' 
which takes any u G Wq ,2 (VL) to the functional A m [it] defined by 

A 2s [u][ip\ = [ A s uA s <pdx, Wtp G Wq' 2 (Q), (2.1) 
Jn 

if m = 2s and 

A 2s+1 [u]{p] = - [ V(A s u) ■ V{A s cp)dx, V<p G iy n ' 2 (O), (2.2) 
Jn 

if m = 2s + 1, where s G No- Thus, the weak formulation of the classic problem 
flUD reads 

(-A)>]M=A(-AHti]M, Vv,Giy n ' 2 (O). (2.3) 

By the Poincare inequality, it follows that the quadratic form defined by 
(— A) n [w][M] for all u G Wq ,2 {VL) is coercive in Wq' 2 (Q), hence the operator (— A) n 
is a linear homeomorphism from Wq ,2 {VL) onto (W n ' 2 (O)) / . Thus the equation 
( 12. 3 p is equivalent to the equation (— A)~ n o (— A) m [tt] = A _1 m, where (—A) - " 
denotes the inverse of (— A) n . It is convenient to endow the space Wq' 2 {Q) with 
the scalar product defined by 

<u 1 ,u 2 > n =(-A) n [u 1 ][u 2 ], (2.4) 
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for all Mi, 1*2 £ Wq' 2 (£1). The norm induced by this scalar product is equivalent 
to the standard Sobolev norm. In the sequel, unless otherwise indicated, we shall 
think of W 

Hilbert space equipped with the scalar product (12 .4p . This 
allows to give a straightforward proof of the following 

Lemma 2.5 Let Q be a bounded open set in M, N and m, n £ N with < m < 
n. The operator Sq = (— A) _n o (— A) m is a non-negative self-adjoint compact 
operator in the Hilbert space Wq' 2 (Q). The spectrum of Sq is discrete and consists 
of a decreasing sequence of positive eigenvalues of finite multiplicity converging to 
zero. Moreover, the equation Squ = /iu is satisfied for some u £ Wq ,2 {Q), /x > 
if and only if equation A2.3\) is satisfied with A = yT l . 

Proof. The equality < SqUi,u 2 > n = (— ^) m [ u i} [^2], for all Ui,u 2 £ Wq' 2 (Q) 
and the symmetry of the operator (— A) m implies that Sq is a self-adjoint oper- 
ator. Since Q is bounded and m < n, the space Wq' 2 (Q) is compactly embedded 
into W™ ,2 {£1). This implies that the operator (— A) m is a compact operator from 
the space Wq' 2 (Q) to its dual. The rest of the proof is trivial. □ 

By Lemma 12.51 and standard spectral theory we deduce the following 

Corollary 2.6 Let Q be a bounded open set in WL N and m, n £ N with < m < n. 
The eigenvalues of problem A2.3\) are positive, have finite multiplicity and can 
be represented as a non- decreasing divergent sequence j £ N where each 

eigenvalue is repeated according to its multiplicity. Moreover, 

\Afl] = A"' m [p] = min max R nm [u], (2.7) 

dimi?=j ' 

for all j £ N ; where R nm [u] is the Rayleigh quotient defined by 

f a \A r u\ 2 dx 

f n \A r u\ 2 dx 
f n \\7A s u\ 2 dx> 

f n \A s u\ 2 dx ' 

L \VA-u\ 2 dx 
J n \VA s u\ 2 dx> 

Clearly, the eigenvalues A™' m [f2] depend on n,m. However, for the sake of 
simplicity, we shall write Xj [O] instead of A"'' m [Q] . 



Rn.m. U 



if n 


= 2r, 


m 


= 28, 


if n 


= 2r, 


m 


= 28 + 1, 


if n 


= 2r + 1, 


m 


= 28, 


if n 


= 2r + 1, 


m 


= 2s + l. 



3 Analyticity results 

Let Q be a bounded open set in R N of class C 1 . In the sequel, we shall consider 
problem (12. 3p in a family of open sets parameterized by suitable diffeomorphisms 
4> defined on Q. Namely, we set 

A« n = \teCm^ N )-- mf 1^-^)1 >q\, 
[ x lt x 2 en \xi — X2\ J 

Xl^X 2 
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where C£(Q;R N ) denotes the space of all functions from O to M. N of class 
C n , with bounded derivatives up to order n. Note that if G Aq then 
is injective, Lipschitz continuous and inf^ |detV0| > 0. Moreover, 0(0) is a 
bounded open set and the inverse map 0*" 1 ) belongs to A^r n y Thus it is nat- 
ural to consider problem (12. 3 p on 0(0) and study the dependence of Aj[0(O)] 
on G Aq. To do so, we endow the space C£(Q;M. N ) with its usual norm 
ll/llq»(n ; R") = sup H < nj xen \D a f(x)\. Note that A% is an open set in C£(Q;R N ), 
see [TTt Lemma 3.11] . Thus, it makes sense to study differentiability and analyt- 
icity properties of the maps i— >■ Aj[0(O)] defined for G Aq. For simplicity, we 
write Aj[0] instead of Aj[0(O)]. As in [17], we fix a finite set of indexes F C N 
and we consider those maps G .Aq for which the eigenvalues with indexes in F 
do not coincide with eigenvalues with indexes not in F; namely we set 

= {0 G Al : A,[0] + A,[0], Vj G F, I G N\F}. 

It is also convenient to consider those maps G A™/^ such that all the eigenvalues 
with index in F coincide and set 

0$? = {0 G ^ : A,J0] = A, 2 [0], V j x ,j 2 G F) . 

For G Ap 1 ^, the elementary symmetric functions of the eigenvalues with 
index in F are defined by 

Af,/J0]= Yl X n[ ( t>]---\M h=l,...,\F\. (3.1) 

h<—<3h 

The main result of this section is the following generalization to poly-harmonic 
operators on smooth domains of the results in [TT| §3] concerning the Dirichlet 
Laplacian on rough domains. 

Theorem 3.2 Let Q be a bounded open set in M> N of class C 1 , n, m G N with 
< m < n, and F be a finite set in N. The set A n /^ is open in C£(Q ; M> N ) and 
the real-valued maps which takes G A n /^ to Ap y h[4>] are real- analytic on A n /^ 
for all h = 1, . . . , \F\. Moreover, if G is such that the eigenvalues \j[</>] 

assume the common value A^[0] for all j G F, and 0(0) is of class C 2n then the 
Frechet differential of the map Ap,h o,t the point is delivered by the formula 

<wm*i = -m p 1 .- 1 ) e ] m) (^) 2 «< ° * < - i) ) • (3.3) 

for allip G C^(0; IR^), where {vi}i(zF is an orthonormal basis in Wq' (0(0)) (with 
respect to the scalar product {2.J$ ) of the eigenspace associated with Ai?[0] ; and v 
denotes the unit outer normal to (90(0). 

Note that formula (13. 3p is a generalization of the celebrated Hadamard for- 
mula, see Grinfeld [11] for a recent paper on this topic; see also Ortega and 
Zuazua [22] for the analysis of associated bifurcation phenomena concerning mul- 
tiple eigenvalues of the biharmonic operator subject to Dirichlet boundary con- 
ditions. 
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In order to prove Theorem 13.21 we consider equation (12.31) on 0(0) and pull it 
back to O. Namely, we consider the equation 

(-A)>]M = A(-A)>M V^e^ o "' 2 (0(O)), (3.4) 

in the unknowns v G W^' 2 (0(O)), A e]0,oo[. Recall that the pull-back to O of 
the classic Laplace operator on 0(0) is defined by 

A^u = (A(u o 0^ -1 ^)) o (3.5) 

for all w G W£ c (0), G The operator A is in fact the Laplace-Beltrami 
operator associated with the change of variables defined by 0. Note that 

A%u = (A s (mo0^ 1 ))) o0 (3.6) 

for all u G W^'^O), G *4|f. For any < m < n, the operator AT can 
be considered as the operator acting from M^' 2 (0) to its dual, which takes any 
u G W-^'^O) to the functional A^[u] defined by 

A;MM = A m [no0(- 1 )][^o0(- 1 )], 

for all <p G W^'^O). More precisely, if m = 2s, s 6 No then 

Al s [u][ v ] = [ AJuAJ^detV^dz, (3.7) 
for all p G W "' 2 (O). If m = 2s + 1, s G N then 

- A}' +1 MM = / V(AJ«) (V0)- 1 (V0)-*V t (A^)|detV0|da;, (3.8) 
Jn 

for all (p G H^ 1 ' (O), where (V0) _1 denotes the inverse of the Jacobian matrix 
of and (V0)~* the transpose of (V0) _1 . Note that the map from W^' 2 (0) to 
W o n ''(0(O)) which maps titoiio 0(- 1 » for all u G W n ' 2 (O) is a linear homeomor- 
phism. Hence, equation (|3.4p is equivalent to 

(-A*)>][d = A(-A,)>M V<^>G<' 2 (0) (3.9) 

where u = v o 0. It is also natural to pull-back the scalar product of Wq' 2 (<J)(Q)) 
to O by setting 

< Mi, U 2 >n,<P = < Ui O (_1) , M 2 ° (_1) >n (3.10) 

for all «i,«2 £ W / o i ' 2 (^)? where < ■,■ > n is the scalar product in W^' 2 (0(O)) 
defined by (12~I]) . By W "£(fi) we denote the Hilbert space PVj 1 ' (O) endowed with 
the scalar product < •, • > n ,<t>- It turns out that the operator S^n) defined in 
Lemma [2751 is unitarily equivalent to the operator defined on Wq'2 (O) by 

T = (-A )-"o(-A^) m . (3.11) 

Thus we can prove the following lemma where C(Wq' 2 (Q)) denotes the space 
of linear bounded operators from W^' 2 (0) to itself and and B s (Wq' 2 (Q)) denotes 
the space of bilinear forms on W^' 2 (0) (both spaces are equipped with their usual 
norms) . 
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Lemma 3.12 Let Q be a bounded open set in M. N of class C 1 , n,m £ No, 
< m < n. The operator defined in A3.11\) is non-negative self-adjoint and 
compact on the Hilbert space Wq^(Q). The equation flff.^p satisfied for some 

v £ Wq' 2 (4>(Q)) if and only if the equation T^u = \iu is satisfied with u = v o 
and fi = A -1 . Moreover, the map from Al to C{Wq'\Q)) x B s {W^\n)) which 
takes £ An to (T^, < •, ■ > n ,<t>) is real- analytic. 

Proof. Since the operator is unitarily equivalent to the operator S^q), 
the first part of the lemma immediately follows by Lemma 12.51 In order to prove 
the real-analytic dependence of upon 0, we note that by standard calculus 

v^-\ / d 2 u du da ri \ 

A 4> u = /, \ a — a n a si + tt- a — a si (3.13) 



for all u £ W 2,2 (fl), where a = (V0) 1 (see also [32J Proposition 3.1]). By for- 
mula f l3TT3|) . it follows that the map from A% x W n ' 2 (fl) to W n - 2 ' 2 {Vt) which takes 
(0, m) £ Aq to A^u is real- analytic. Thus also the maps from Aq x W n ' 2 (Q) to 
L 2 {9) which take (0, u) £ Aq to A^u for all s £ N with < s < n/2, are 
real-analytic since they are compositions of real-analytic maps. This, combined 
with formulas (13. 7p and (13. 5J) , implies the real-analytic dependence of and 
< •, • > n ^ upon 0. □ 

Proof of Theorem 13.21 We denote by Hj[4>\, j £ N, the eigenvalues of T^. 
By Lemma 13.121 fJ,j[(f>) = Aj 1 [0] for all j £ N, hence the set ApJ^ coincides with 
the set {<f>eA%: Hj[<j>\ ± m[<j)\, Vj £ F, / £ N \ F}. By Lemma GD2J T^ is self- 
adjoint with respect to the scalar product < -, • > n ,<f> an d both and < *, • > n ,<i> 
depend real-analytically on 0. Thus, by applying [T71 Thm. 2.30], it follows that 
ApQ is an open set in C£(Q ; WL N ) and the functions which take £ Ap^ to 

r F M= E »M---»iM (3-i4) 

ji,...,j h eF 
j .a 

are real-analytic for all h = 1, . . . , \F\. Since 

A ^[0J = ~F FTT' ( 3 - 15 ) 



r 



for all h — 1,..., where r^ jO [0] = 1, it follows that k F ^ h [(f>] depends real- 



analytically on £ A^q. 



It remains to prove formula (13. 3p . Let £ O^'q, \f[4>] an d {^}«sf be as in 
the statement. We set ui = v \ o for all / £ F and we note that {u{\i & f is an 



rn,2, 

o,4> \ 

A^ x [0] of the operator TV. By [I?] Thm. 2.30], it follows that 



orthonormal basis in W™'~ (f2) for the eigenspace corresponding to the eigenvalue 



dl^iVM = x F~ h [4>]( {F h l _ S) E < d l=^MN.«i >m (3-16) 

^ ' leF 

for all -0 £ (7^(12 jR^). Note that by standard regularity theory (see e.g., Ag- 
mon PQ Thm. 9.8]) v l £ VT 2n ' 2 (0(fi)) for all I £ F. 
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By standard calculus, equalities (12 .4p . (I3.10p . Theorem 15.71 by observing that 
} = on d0(O) and (-A) n ^ = \ F [<p}(-A) m v h we have 



= (d|^(-A^) m M)NM - A^[0](d|^(-A )"[^])Nh] 
= - / (f^r) C-uda-2 [ (-A)"^ ■ (da 

I C-vda + 2A F 1 [0] / (-A)"vjVv, • 

2 



A 7 /[o] / (^) (3/1.7) 



9</>(f}) 



where we have set C — "0 ° ^- Formula (13. 3D easily follows by combining for- 
mulas fl3TT5l -f l3T7l) . □ 



4 Isovolumetric perturbations 

Consider the following extremum problems for the symmetric functions of the 
eigenvalues 

min A Fth [4>] or max A F>h [(f)], (4.1) 

V[4>]=ccmst V[4>]=ccmst 

where V[0] denotes the iV-dimensional Lebesgue measure of 0(0). 

Note that if G Aq is a minimizer or maximizer in (14. ID then is a critical 
domain transformation for the map i— >■ Ap^[0] subject to volume constraint, 
i.e., 

Ker d|^V C Ker d^A^, (4.2) 

where V is the real valued function defined on Aq which takes G A^ to V[(p]. 

The following theorem provides a characterization of all critical domain trans- 
formations 0. See [TB] for the case of the Dirichlet and Neumann Laplacians. 

Theorem 4.3 Let Q be a bounded open set in M> N of class C 1 , n, m G N with 
< m < n, and F be a finite subset of N. Assume that G B^'n such 
that 0(0) o/ c/ass C 2n and i/iai £/ie eigenvalues Aj[0] have the common value 
\ F [4>] for all j G F. Let {vi}i eF be an orthornormal basis in Wq' 2 (0(0)) of the 
eigenspace corresponding to X F [(f>}. Then (ft is a critical domain transformation 
for any of the functions A F ^, h = 1, . . . , \F\, with volume constraint if and only 
if there exists Cel such that 

2 

= C, on 90(0). (4.4) 

Proof. Note that V[<f>] = J n |detV0|dx, hence by formula (I5.10p it follows 
that 

d\t=iV[1>)= f (J, o . vda, (4.5) 
Jd4>(n) 
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d n V[ 



dv r 



for all xjj G C^ip, ; M. N ). The proof of (I4.4p follows immediately by formulas ( 13. 3 p 
and (14. 5p . and by observing that condition (14.21) is satisfied if and only if there 
exists c G R (a Lagrange multiplier) such that 

d|^A^ = cd|^y. (4.6) 

□ 

Finally, we can prove the following 

Theorem 4.7 Let the same assumptions of Theorem \4-3\ hold. If <f)(Q) is a ball 
then condition ^4-4\) i> s satisfied. 

Proof. Without any loss of generality, we assume that 4>(fl) is a ball with 
radius R centered at zero. By the rotation invariance of the Laplace operator, 
{vioA}i & f is an orthonormal basis of the eigenspace corresponding to Xf[4>] f° r all 
A G Otv(R) where On(R) denotes the group of orthogonal linear transformations 
in R^. Since both {vi}i & p and {vi o A}i e p are orthonormal bases of the same 

space, it follows that Xj=i vf o A = Xj=i v f > f° r aU A G O n (R). Thus Ya=i v l * s 
a radial function. By differentiating with respect to the radial coordinate r, by 
the Leibniz formula and by recalling that all derivatives up to order n — 1 of the 
eigenfunctions vanish at the boundary of <f>(£Y), we obtain that 



dr 



2n 



d 2n vf ^ {2n\ fd k vid 2n k vi . 



r=R k=0 



k 7 V <9r fc <9r 



2n-fc 



2n\ / <9""i^ 



n I \ dr 



r=R 



Since X]/gf dr 2 " ^ s a ra( iial function, then by formula ( 14. 8ft we conclude that the 



EieF&) 2 is constant on □ 



It would be interesting to clarify whether balls are local minimizers or max- 
imizers for the eigenvalues or their symmetric functions. With regard to this, 
we mention that it is proved in Wolf and Keller [24"| Thm. 8.3] that the circular 
disk in the plane is a local minimizer for the third eigenvalue of the Dirichlet 
Laplacian. 

5 A formula for the Frechet differential of the 
'poly-Laplace-Beltrami' operator 

In this section we prove Theorem 15 . 71 which has its own interest since it provides 
an explicit formula for the Frechet differential with respect to of the weak 
'poly-Laplace-Beltrami' operator A 1 ^ defined by ( 13. 7p . ( 13. 8p . That formula has 
been used in the proof of (13. 3p . 

Lemma 5.1 Let Q be a bounded open set in M, N of class C 1 , s G N, U\ G L 2 (Q), 
U2 G Wq S,2 (Q). Let <j) G Aq and Vi = Ui o <p^ l \ i = 1,2. Assume that </>(0) is of 
class C 1 and that v x G W 2s ' 2 (4>({l)), v 2 G W 2s+1 ' 2 (4>({l)). Then 

[ u a d|^AJ[VHdetV0|dz 
Jn 

= I (^iVA s t; 2 - A s Vl Vv 2 ) -£dy- [ v 1 A s v 2 ( ■ uda, (5.2) 
for all if) G Cl s {n ; R N ), where ( = ip o 
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Proof. First, we recall the following formula from [T%| Lemma 3.42] which 
holds for any u G W 2,2 (Q): 

(d|^[V»]«) o <f>™ = -2j2 it ~ E 5? A • (5-3) 



We observe that 



s-1 



E ^ 



(5.4) 



h,fc=0 



By formulas (15.31) and ( 15.4ft . by changing variables in integrals and integrating 
by parts, we obtain 



uid|^Ai[^]u 2 |detV0|da; 



/l+fe=S— 1 
5-1 



E L. E 



f x 9y» ^ % / 



<9A% dA k v 2 fdQ d(, 



h,k=0 - / ^(°) i,j=l 

h+k=s-l 

h n , \k+l„, i r?A/i.. V7 A fc„ 



+ 



-(A%A fc+1 ?; 2 + VA%!VA%)divC dj/, (5.5) 
see also [HI Formula (3.45)]. Moreover, integrating by parts yields 



+ / ASVA fc+1 va • Cdy + / A^iA fc+1 w 2 divC^ 

+ / VA% • VA fc w 2 divCc/?/ - / A ft+ %VA% • Cdn. (5.6) 
j^(n) «/<£(«) 

By observing that the first summand in the right-hand side of (15. 6p vanish if 
k < s — 1, and by combining ( 15. 51) and ( 15. 61) . we obtain a telescopic sum and we 
deduce the validity of ( 15.21) . □ 

Theorem 5.7 Let Q be a bounded open set in W N of class C 1 , n G N, U\,v,2 G 
W / ^' 2 (n). Lei G and t> j = Mj o c^" 1 ), i = 1,2. Assume that 4>{Sl) is of class 
C 1 and that v u v 2 e W 2n ' 2 (4>(tt)) . Then 

((-A) n Vl Vv 2 + (-A) n i; 2 Vt;i) ■ Cdy,(5.8) 
/or a// V e C b n ( fi ! rN )> w/iere C = ip ° </> M) ■ 
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Proof. First, we consider the case where n is an even number of the form 
n = 2s with s G No. By formula (13. 7p and standard calculus we have 



n 



(d|^Aj'MM) (u 2 ) = I dl^AJti^^AJrial det D~^\dx 



n 



+ ! AJtnAJuadl^ldetD^ltVldx. (5.9) 

J f2 

Moreover, by standard calculus 

(d|^ (detV0) [</>]) o detV^" 1 ) = div U o \ , (5.10) 



hence 



/ Al Ul Alu 2 d\^i\ det D(j>\[ij;}dx= I A^A^divC^. (5.11 
Jn J $(a) 



Formula ( 15. 8 p easily follows by combining formulas ( 15.2ft , ( 15. 9 ft , ( 15. lift , by integrat 
ing by parts and by observing that A s Vi = |^ on 90(0) since t>j G Wq S ' 2 (^(^)) 
Now, we consider the case where n is an odd number of the form n = 2s + '. 
with s e N . By formula (13. 8 p and standard calculus we have 



{d\^Af +1 m Ul )) (u 2 ) = f VA S ^(VC + W'QV'A^dy 

J4>(U) 

f VA s Vl VA s v 2 dw(dy - I VA s Vl V Ud\^%u^\) o dy 

- f V ((dl^AJu!^]) o 0^)) VA^ 2 dy.(5.12) 



Moreover, integrating by parts yields 

/ VA S ^(VC + V t ()V t A s v 2 dy 

Ji(Q) 



N 



I / dCh dA s Vl dA s v 2 | d(k dA s v 1 dA s v 2 \ ^ 

^— ' /t,„, \ r)-?/, r)-?/, f)n, rhi, fin, fin, I 



h,k=l 



\9yk 9y h dy k dy h dy h dy k 
dA s vi dA s v 2 



90(0) dv dv 



■ vda 



A r (d 2 A s v 1 dA s v 2 ^ + dA s Vl d 2 A s v 2 ^ ^ 
2-^ l~<~\dy h dy k dy k h dy h y\ h l 



h,k=l 
N 



0(C) 



, /d 2 A s <uidA s <u 2 ^ + dA s v l d 2 A s v 2 ^ j 



h,k=l 



(>A '' { <>A U2 (-vda+ I VA s VxVA s v 2 divCdy 



90(n) <9z/ ^ J<j>(n) 



-[ (A s+1 v 1 VA s v 2 + A s+1 v 2 VA s v 1 )-(dy. (5.13) 
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By integrating by parts, changing variables in integrals and using formula 
( 15. 2p . we obtain 

/ VA^v((rf|^A> J M)o0(- 1 ))^ 

A^d|^A^#]|detV0|cfe 
= - / (A S+1 ^VA S ^- - A 2s+1 ^V^) ■ Qdy (5.14) 

for all i,j G {1,2}. Finally, formula (I5.8P easily follows by combining formu- 
las dSH2D-(EIHD and by observing that = on d<j)(£l) since v { e 
W" o 2s+1 ' 2 (0(O)). " "' □ 

Acknowledgments: The second author acknowledges financial support from 
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